EXPERIMENTAL  DETERMINATION  OF 
THE  MOLECULAR  DIFFUSION  COEFFICIENT  OF  GASES 
BY  ENHANCED  DISPERSION  IN  OSCILLATORY  FLOWS 


By 

GREGORY  ALLEN  ZACHRICH 


A DISSERTATION  PRESENTED  TO  THE  GRADUATE  SCHOOL 
OF  THE  UNIVERSITY  OF  FLORIDA  IN  PARTIAL  FULFILLMENT 
OF  THE  REQUIREMENTS  FOR  THE  DEGREE  OF 
DOCTOR  OF  PHILOSOPHY 

UNIVERSITY  OF  FLORIDA 


1995 


ACKNOWLEDGEMENTS 


Dr.  Ulrich  H.  Kurzweg,  chair  of  my  committee,  lead  me 
through  this  research  project.  Although  he  performed  many 
required  tasks,  specifically  I thank  him  for  reviewing  this 
paper's  various  preliminary  versions  suggesting  several 
modifications  to  improve  the  final  copy.  Also,  his 
continued  effort  to  keep  me  focused  on  the  main  topic  of 
this  dissertation  is  much  appreciated. 

Dr.  E.  Rune  Lindgren,  whom  I consider  a friend, 
prompted  me  to  further  my  education.  I thank  him  for 
allowing  me  to  discuss  many  of  my  thoughts  in  confidence 
with  him. 

Dr.  Marc  Jaeger  was  extremely  helpful  in  providing  his 
laboratory  equipment  and  assisting  in  the  collection  of 
data.  Without  his  cooperation  this  work  would  not  have  been 
possible . 

I thank  the  remainder  of  my  committee,  Drs.  Shyy  and 
Lear,  for  serving  in  their  respective  capacities  as  member 
and  external  member. 

I truly  appreciate  the  valued  advice  and  guidance  of 
Shirley  Robinson,  who,  among  other  things,  reminded  me  of 
the  numerous  deadlines  that  arose. 


11 


My  mother,  father,  brothers  and  sister  provided 
unquant if iable  moral  support  during  this  long  endeavor,  for 
which  I am  most  grateful. 


iii 


TABLE  OF  CONTENTS 

page 

ACKNOWLEDGEMENTS ii 

LIST  OF  TABLES vi 

LIST  OF  FIGURES vii 

LIST  OF  SYMBOLS ix 

ABSTRACT xii 

CHAPTER  1 INTRODUCTION 1 

CHAPTER  2 ESTIMATING  THE  DIFFUSION  AND  VISCOSITY 

COEFFICIENT  IN  GASES 5 

Estimating  the  Diffusion  Coefficient 5 

Estimating  the  Viscosity  Coefficient 9 

CHAPTER  3 ANALYTICAL  SOLUTION  FOR  ENHANCED  DISPERSION 

IN  OSCILLATING  FLOW 14 

Velocity  Profile  Solution 14 

Concentration  Solution 23 

Determining  the  Enhanced  Diffusion  Coefficient. . . 26 

Tidal  Displacement  Solution 32 

Enhanced  Diffusion  Coefficient  in  Final  Form 38 

Graphical  Presentation 39 

CHAPTER  4 EXPERIMENTAL  PROCEDURE 41 

Enhanced  Diffusion  Process 41 

Experimental  Procedure 44 

Experimental  Assumptions 48 

Diffusion  Coefficient  is  Independent  of 

Concentration 48 

Incompressible  Gases 49 

Linear  Concentration  Gradient 50 

Pressure  Gradient  is  a Function  of  a Pure 

Sine  Wave  in  Time 50 

Cross-Sectional  Velocity  Profile  is 

Fully  Developed 53 

Transverse  Diffusion  Effects  Overwhelms 

Axial  Diffusion 53 

Flow  is  Completely  Laminar 54 

Steady  State  Conditions  have  been  Achieved. . 54 

Tube  is  Completely  Straight 55 


IV 


All  Measurements  are  Accurate 55 

No  Leaks  To  or  From  Tube  System 56 

Tidal  Displacement  Remains  Small 57 

Deriving  the  Non-dimensional  Expressions  Used  in 
the  Experimental  Measurements  of  the  Two  Setups..  57 
Determining  Density  and  Viscosity  Coefficients 
Experimentally 63 

CHAPTER  5 EXPERIMENTAL  RESULTS 64 

Setup  One 64 

Setup  Two 69 

Pressure-Displacement 73 

CHAPTER  6 CONCLUSIONS 75 

Concluding  Remarks 75 

Future  Work 77 

APPENDIX 79 

TABLE  OF  REFERENCES 84 

BIOGRAPHICAL  SKETCH 8 7 


v 


LIST  OF  TABLES 


Table  page 

2.1  Tabulated  Values  of  Parameters  for 

Various  Gases 8 

2.2  Predicted  and  Reference  Diffusion  Coefficient  8 

2.3  Predicted  and  Reference  Viscosity 11 

2.4  Diffusion  Collision  Integral 12 

2.5  Viscosity  Collision  Integral 13 

5.1  Raw  Data  Setup  1 67 

5.2  Results  Setup  1 68 

5.3  Raw  Data  Setup  2 70 

5.4  Results  Setup  2 71 

5.5  Pressure-Displacement  Results 74 


vi 


TABLE  OF  FIGURES 

Figure  page 

3.1  Plot  of  Log  J0(  (-i)  5 x) 19 

3.2  Plot  of  Log  J,((-i)  5 x) 20 

3.3  Plot  of  ber0,  bei0,  ber, , and  bei, 21 

3.4  Velocity  Profiles  at  0 Degree 22 

3.5  Velocity  Profiles  at  30  Degree 22 

3.6  Velocity  Profiles  at  60  Degree 22 

3.7  Velocity  Profiles  at  90  Degree 22 

3.8  Velocity  Profiles  at  120  Degree 22 

3.9  Velocity  Profiles  at  150  Degree 22 

3.10  Pressure-Displacement  Relationship 

(Smaller  a) 35 

3.11  Pressure-Displacement  Relationship 

(Larger  a) 36 

3.12  Phase  Between  Pressure  Gradient  and  Velocity  37 

3.13  Contour  Plot  of  X = Deff  /w  Ax2  verses  a and  S 40 

4.1  Schematic  of  Setup  1 46 

4.2  Schematic  of  Setup  2 46 

4.3  Photograph  of  Setup  2 47 

4 . 4 Photograph  of  Linear  Motor  Used  to  Produce 

Sinusoidal  Axial  Displacements  of  Fluid 
within  Conduit 50 

4.5  Photograph  of  Pressure  Variation  and 

Tidal  Volume 51 


vii 


4.6  Contour  Plot  of  N(a,S) 61 

5.1  Experimental  Data  Setup  1 69 

5.2  Experimental  Data  Setup  2 71 

5.3  Experimental  Data  Pressure-Displacement 

Relationship 74 


vm 


KEY  TO  SYMBOLS 


a Womersley  Number,  a = a*(oj  / v)  5 (non-dimensional) 
7 Concentration  Gradient  (g/cm4) 

Ax  Tidal  Displacement  (cm) 

e Intermolecular  Distance  (cm'8  = A) 

r / Radial  Distance,  rj  = r/a  (non-dimensional) 

6 Angular  Direction  (radian) 

X Dispersion  Function  (non-dimensional) 

M Dynamic  Viscosity  (g/cm  s) 

/iAIR  Dynamic  Viscosity  of  Air  (1  cP  = 0.01  g/cm  s) 
v Kinematic  Viscosity  (cm2/s) 

p Density  (g/cm3) 

Atomic  Diffusion  Volume  Increments 
aAB  Intermolecular  Distance  (A) 

a;  Angular  Frequency  (1/s) 

Collision  Integral  for  Viscosity 
^d,ab  Collision  Integral  for  Diffusion  Coefficient 

a Radius  of  Tube  (cm) 

A Cross-Sectional  Area  (cm2) 

c Acoustic  Velocity  (cm/s) 

C Concentration  (g/cm3  or  %) 

Dab  Molecular  Diffusion  Coefficient  (cm2/s) 


IX 


Deff  Effective  Enhanced  Diffusion  (cm2/s) 

Dm  Molecular  Diffusion  Coefficient  (cm2/s) 
dS  Circumf erencial  Distance  dS  = rj  dd  (non-dimensional) 
f(r?)  Non-dimensional  Velocity  Function 
F(rj)  Non-dimensional  Concentration  Function 
g(r?)  Non-dimensional  Concentration  Function 
G(rj)  Non-dimensional  Concentration  Function 
h (rj)  Non-dimensional  Concentration  Function 
J0  Bessel  Function  of  First  Kind  of  Zero  Order 
Bessel  Function  of  First  Kind  of  First  Order 
kb  Boltzman's  Constant  (1.38054  x 10"23  J/molecule  °K) 

K Magnitude  of  Pressure  Grad.  / density  (Atm  cm2/g) 

L Length  of  Tube  (cm) 

M Molecular  weight  (g/gmole) 

N Non-dimensional  Function 

P Pressure  (Atmospheres) 

Q Volumetric  Flow  Rate  (cm3/s) 
r Radial  Direction  (cm) 

R Universal  Gas  Constant  (8.3143  J/gmole  °K) 

S Schmidt  Number,  S = v / DM  (non-dimensional) 
t Time  (s) 

T Temperature  (°K) 

Tc  Critical  Temperature  (°K) 

u Velocity  (cm/s) 

U0  Mean  Velocity  (cm/s) 

V Volume  of  Isolation  Chamber  (cm3) 


x 


Vc  Critical  Specific  Volume  (m3/Kg) 
x Axial  Direction  (cm) 

Y0  Bessel  Function  of  Second  Kind  of  Zero  Order 


xi 


Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Doctor  of  Philosophy 

EXPERIMENTAL  DETERMINATION  OF 
THE  MOLECULAR  DIFFUSION  COEFFICIENT  OF  GASES 
BY  ENHANCED  DISPERSION  IN  OSCILLATORY  FLOW 

by 

Gregory  Allen  Zachrich 
May  1995 

Chairman:  Ulrich  H.  Kurzweg 

Major  Department:  Aerospace  Engineering, 

Mechanics,  and 
Engineering  Sciences 

This  paper  presents  an  experimental  study  of  an 
alternative  method  to  determine  the  molecular  diffusion 
coefficient  utilizing  an  enhanced  axial  dispersion  process. 
Enhanced  dispersion  is  a process  by  which  a species  in  a 
fluid  oscillating  in  a tube  can  experience  a vastly 
amplified  rate  of  effective  diffusion  as  compared  to  normal 
molecular  diffusion.  This  enhancement  is  promoted  by  an 
interaction  between  the  boundary  layer  and  the  core  within 
conduits  as  the  fluid  oscillates.  The  species  is 
transferred  first  by  diffusing  radially  into  the  boundary 
layer  and  later  in  the  cycle  diffusing,  again  radially,  back 
into  the  oscillating  core  of  the  fluid,  being  transported 
axially  a large  distance  from  the  originating  location. 
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This  enhanced  dispersion  process  is  completely 
described  mathematically  by  a function  dependent  on  the 
properties  of  the  kinematic  viscosity  of  the  carrier  fluid, 
the  molecular  diffusion  coefficient  of  the  diffusing  species 
into  the  carrier  fluid,  the  tube  radius,  the  angular 
freguency  of  oscillation,  and  the  axial  tidal  or  bulk 
displacement . 

Measurements  of  the  diffusional  flux  were  obtained  with 
a sinusoidal  freguency  varying  between  1.0  and  10.0  Hz, 
tidal  displacement  extending  from  13.2  to  50.9  cm,  and  tube 
radius  ranging  from  0.19  to  1.03  cm  in  time-periodic  steady 
state  conditions.  The  experiment  was  conducted  on  two 
separate  setups,  one  a system  of  multiple  capillary  tubes 
connecting  two  reservoirs  with  no  net  flow  in  the  tubes,  the 
other  consisting  of  a single,  larger  tube  between  a 
reservoir  and  the  atmosphere  with  a net  flow.  All 
measurements  were  made  at  room  temperature  (-300°  K)  and  at 
one  atmosphere  pressure.  Nine  experimental  runs  were 
completed  with  nitrogen  diffusing  into  oxygen  and  seven  of 
oxygen  diffusing  into  carbon  dioxide.  The  experimental 
molecular  diffusion  coefficient  was  determined  by  a least 
sum-of-squares  of  error  between  the  theoretical  curve  and 
the  collected  data.  Also  a predicted  'tuning  effect'  in 
this  dispersion  process  was  verified  with  optimum  effective 
diffusion  occurring  in  02-C02  mixtures  near  the  Womersley 
number  of  4 . 
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The  molecular  diffusion  coefficient  of  oxygen  into 
carbon  dioxide  was  determined  to  be  0.136  cm2/s  when  reduced 
to  0°  C and  one  atmosphere  of  pressure,  compared  to  a 
tabulated  value  of  0.139  cm2/s  at  the  same  conditions. 
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CHAPTER  1 


INTRODUCTION 

Diffusion  is  a process  of  mass  transport  due  to  the 
random  molecular  motion  of  the  individual  mass  particles. 
Although  Graham  made  experimental  studies  of  the  phenomena 
as  early  as  1803,  the  mathematical  study  of  diffusion  was 
first  investigated  by  Fick  in  1855.  He  found  that  diffusion 
was  governed  by  the  same  mathematical  expressions  as  Fourier 
had  derived  for  heat  conduction  in  1822.  Therefore,  many 
solutions  found  for  heat  conduction  are  directly  applicable 
to  the  diffusion  processes  (Crank  1975) . 

The  study  of  diffusion  in  gases  was  first  recorded  in 
1870  when  Loschmidt  devised  a method  to  experimentally 
determine  the  diffusion  coefficient.  His  apparatus 
consisted  of  two  cylinders  placed  vertically,  separated  by  a 
stopcock.  Into  this  device  were  placed  the  two  gases  that 
were  to  be  studied.  The  more  dense  gas  was  placed  in  the 
lower  chamber  to  preclude  the  likelihood  of  convective 
currents  internally,  causing  mixing  that  would  artificially, 
but  erroneously,  increase  the  determined  numeric  diffusion 
coefficient.  Once  the  experiment  was  prepared,  the  stopcock 
was  opened,  and  after  a predetermined  time  was  closed. 

After  the  experiment  was  concluded,  the  chemical  composition 
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of  the  two  chambers  were  determined,  usually  by  chemical 
adsorption  of  one  of  the  gases,  which  provided  the  raw  data 
to  determine  the  diffusion  coefficient. 

In  1907  an  improvement  to  the  original  method,  a 
sliding  disk  that  had  the  full  cylinder  diameter  available 
for  the  diffusion  was  implemented.  This  allowed  for  a 
greater  rate  of  diffusion,  due  to  the  increased  area,  which 
increased  the  accuracy  of  the  measurement.  In  this 
apparatus,  temperature,  pressure  and  concentration 
dependance  were  investigated. 

It  was  determined  that  the  temperature  dependance  went 
at  T1 5 power.  The  pressure  dependance  was  found  to  be 
inversely  proportional  to  the  pressure,  and  the 
concentration  dependance  was  generally  very  limited,  but 
nearly  the  square  root  of  the  molecular  weight  of  the  gas 
mixture. 

The  study  of  accelerated  diffusion  in  pipe  flow  was 
investigated  by  Taylor  (1958),  who  studied  both  laminar  and 
turbulent  steady  flow  conditions.  As  an  experimental 
verification  of  his  theory,  he  injected  salt  into  flowing 
water  and  obtained  a diffusion  distribution  some  distance 
downstream.  This  study  invalidated  earlier  assumptions  that 
the  salt  would  flow  as  a 'block'  downstream,  and  proved  that 
there  was  an  interaction  between  the  boundary  layer  and  the 
core  flow.  Taylor  (1953)  showed  mathematically  that  the 
apparent  dispersion  coefficient  could  be  written  as 
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U 2 ,3  2 

K=  , where  U is  the  mean  flow  velocity,  a is  the  tube 
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radius,  and  DM  is  the  molecular  diffusion  coefficient. 

Aris  in  1956  showed  that  Taylor,  by  neglecting  the 
axial  diffusion  in  the  derivation  of  his  equation,  had  to 
place  many  restrictions  on  the  applicability  of  his  result. 

A more  precise  statement  is  that  the  axial  molecular 

• • 77 2 r)  2 

diffusion  rate  remained  in  the  equation  as  K = DM+ . 

"48  Dm 

This  eliminated  all  of  the  restrictions  of  Taylor  (Aris, 
1956)  . 

Taylor  also  showed  that  the  molecular  diffusion 
coefficient  can  be  determined  by  measuring  the  spread  caused 
by  this  enhanced  diffusion  method  (Taylor,  1953) . Since 
that  time,  the  measurement  of  the  diffusion  coefficient  has 
been  improved  through  a method  of  streaming  diffusion,  which 
is  essentially  identical  to  the  method  described  by  Taylor. 
This  method  uses  a technique  of  injecting  a small  amount  of 
fluid  into  a laminar  flow  of  a carrier  fluid.  The  smaller 
quantity  is  dispersed  throughout  the  carrier  fluid  by  the 
action  of  radial  diffusion.  This  spread  can  then  be  used  to 
determine  the  molecular  diffusion  of  the  generally  smaller 
quantity  into  the  generally  larger  quantity  carrier  fluid. 
The  advantage  in  this  method  is  the  amplification  factor 
that  is  achieved  by  this  spreading  action.  This  method  is 
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now  generally  called  the  'point  source'  method  for 
determining  molecular  diffusion  coefficient. 

Harris  and  Goren  (1967)  first  examined  axial  dispersion 
in  capillaries  under  pure  oscillatory  conditions.  This  was 
followed  by  Watson  studying  this  oscillatory  diffusion  in 
much  greater  detail  in  1982.  Watson  derived  the  generalized 
equations  for  enhanced  axial  dispersion  for  any  given  tube 
cross-section  and  also  with  superimposed  steady  net  flow. 

Various  papers  have  been  written  on  the  subject  of 
enhanced  diffusion  by  oscillatory  motion  (Dreyer  et  al. 

1968)  (Joshi,  et  al.,  1983)  (Jaeger  and  Kurzweg,  1983) 
(Kurzweg  and  Jaeger,  1986,  1987)  (Jaeger  et  al.  1991) 

(Jaeger  and  Kalle,  1992).  These  papers  generally  proved 
experimentally  that  the  theory  was  in  good  agreement  with 
experimental  data. 

Because  of  the  practical  limitation  of  experimentally 
obtaining  diffusion  coefficients  for  all  conditions  and  all 
possible  combinations  of  components,  a method  was  needed  to 
reliably  estimate  the  diffusion  coefficient.  This  was 
attempted  by  Chen,  Othermer  and  others  starting  in  the  early 
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CHAPTER  2 


ESTIMATING  THE  DIFFUSION  AND  VISCOSITY 
COEFFICIENT  IN  GASES 

The  actual  documented  data  on  molecular  diffusion  is 
quite  limited.  It  is  rather  unlikely  that  the  actual 
conditions  required  for  a specific  request  of  a molecular 
diffusion  coefficient  has  experimentally  been  determined. 
Because  of  this,  much  work  has  been  performed  on  devising  a 
method  to  estimate  the  diffusion  coefficient  for  any  given 
components,  temperature  and  pressure. 

Estimating  the  Diffusion  Coefficient 

Most  of  the  methods  for  estimating  the  diffusion 
coefficient  require  the  use  of  tables  of  self-diffusion. 
These  tables  are  based  on  identical  and  undistinguishable 
components  diffusing  into  one  another.  These  tables  are 
normally  generated  by  using  isotopes  of  the  various  gases  to 
closely  approximate  the  identical  and  undistinguishable 
requirement.  From  these  test  runs,  several  parameters  are 
calculated  such  as  the  collision  diameter,  an  intermolecular 
characteristic  distance,  and  a collision  integral.  Kinetic 
gas  theory  and  quantum  mechanics  provide  these  parameters 
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which  may  subsequently  be  used  to  estimate  the  diffusion  of 
one  substance  into  another  substance. 

Reid  and  Sherwood  (1966)  recommended  using  the 
following  equation  as  an  estimator  of  the  molecular 
diffusion  coefficient. 

n _ 0 .0018583r3/2  I 1 , 1 ,, 

" P^ab)2^d,ab  N Ma 

In  this  equation  T is  the  temperature  in  degrees  Kelvin,  P 
is  the  pressure  in  Atmospheres,  and  a ^ is  the  molecular 

diameter  in  Angstroms,  which  is  determined  from  quantum 
mechanics  and  is  tabulated  for  various  gases.  The  collision 
integral , 0DAB , is  determined  from  a tabulation  (Table  2.4) 

of  the  temperature  forming  a non-dimensional  term  kgT/e^ 

where  kB  is  Boltzman's  constant,  the  Lennard-Jones  force 

constant,  e^,  is  estimated  as  e ^ = JeaeB  and  a ^ is  estimated 

as  °ab  = ^ (°a  + °b>  • Again  eA , eB , aA , and  aB  are  tabulated 

for  the  various  gases.  The  molecular  weights  of  the  gases 
are  MA  and  MB.  For  equation  2.1  the  average  error  for  114 

measured  values  of  diffusion  coefficient  in  65  binary 
systems  at  various  temperatures  was  7.5%  (Skelland  1974). 

Also,  Fuller,  Schettler,  and  Giddings  (1966)  suggested 
using  the  equation 
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o . ooioor7/4  _i_  + _i 


(2.2) 


P[(Ev)i/3+(Ev)y3]2N  ^ 


Here  £v  is  the  sum  of  Vi  which  are  tabulated  Atomic 

Diffusion  Volume  Increments.  This  equation  has  an  average 
error  of  6.9%  for  114  measured  values  of  65  binary  systems 
(Skelland  1974) . 

Chen  and  Othmer  (1962)  suggested  the  version, 


where  P is  in  atmospheres,  Tc  the  critical  temperature  in 

degrees  Kelvin,  Vc  the  critical  specific  volume  in  cubic 

centimeters  per  gram-mole.  This  equation  has  an  average 
error  of  8.5%  for  114  measured  values  of  65  binary  systems 
at  various  temperature  (Skelland  1974) . 

Othmer  and  Chen  (1962)  suggested  the  equation 


(2.3) 


1.23 


(2.4) 


Here  | xalz  is  the  dynamic  viscosity  measured  in  cP 


(CentiPoise) , where  1 cP  = 0 . 01 g /( cm  • s)  , at  the  temperature 
and  pressure  of  the  system  in  question.  This  equation  has 
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an  average  error  of  12.6%  for  the  same  measured  values 
(Skelland  1974) . 

At  room  temperature  (300°  K)  and  one  atmosphere  of 
pressure,  the  tabulated  values  are  shown  in  Table  2.1. 
Table  2.1  Tabulated  Values  of  Parameters  for  Various  Gases 


Gas 

V 

cm2 

s 

P 

kg 

i77 3 

0 

in'10 

e/k 

°K 

Vl 

°K 

cm3/ g 

Air 

0.1566 

1.177 

3.711 

78.6 

20.1 

132.4 

2.857 

C02 

0.085 

1.788 

3.941 

195.2 

26.9 

304.2 

2 . 174 

He 

1.240 

0.1624 

2 . 551 

10.22 

2.88 

5.25 

14 .43 

n2 

0.1553 

1.138 

3.798 

71.4 

17 . 9 

126.1 

3.215 

02 

0.1585 

1.299 

3.467 

106.7 

16.6 

154.4 

2.326 

v,p  (Edwards  et  al .)  o , e/k,  Vi  (Skelland)  Tc,  Vc  (Perry  & Chilton) 


At  3 00 °K  and  one  Atmosphere,  the  following  table 
provides  the  calculated  estimates  for  various  combinations 
of  gases. 

Table  2.2  Predicted  and  Reference  Diffusion  Coefficient 


Gas 

Eq  2.1 

cm2/ s 

Eq  2.2 

cm2/ s 

Eq  2.3 
cm2/ s 

Eq  2.4 

cm2/ s 

Adjusted 

Table 

cm2/ s 

Original 

Table 

cm2/ s 

He->C02 

0 .828 

0.801 

1.063 

1.03 

02->C02 

0.154 

0.163 

0.168 

0.171 

0.164 

0.139 

N2->02 

0.208 

0.210 

0.217 

0.201 

0.213 

0.181 

Reference  values  (Hirschf elder , et  al.  1949) 


The  table  reference  has  been  adjusted  to  300  degree  Kelvin 
from  0 degree  Centigrade  by  the  7/4  power  rule.  (Chapman  & 
Cowling  1953) 
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Estimating  the  Viscosity  Coefficient 

The  dynamic  viscosity  of  a low  density  gas  can  be 
predicted  by  the  equation 

\x  = 2 . 6693.x  10~5  , (2.5) 

o2 

where  |i  is  in  g/{cm-sec)  , M is  the  Molecular  Weight  of  the 
gas,  T is  the  Temperature  in  degrees  Kelvin,  o is  in 
angstroms,  which  is  a tabulated  value,  and  0^  is  also  a 

kT 

tabulated  value  (Table  2.5)  determined  from  as  in  the 

diffusion  coefficient  prediction  equations,  but  of  a 
different  value  (Bird,  et  al.  1960). 

The  National  Bureau  of  Standards  have  produced  formulas 
for  predicting  the  dynamic  viscosity  of  several  gases. 

These  are  simply  interpolation  formulas  from  a best  fit 
quadratic  curve  in  both  temperature  and  pressure  to 
tabulated  data. 

\xaiz  = 17  0.256  + 0.6  05434  T 
- 0 . 001332  r2  + l . 17237  P (2.6) 

+ 0 . 125639  P2 

= 167 . 214  + 0 . 392728  T 

+ . 00122474  T2  + 1 . 24165  P 
+ 0 . 0987206  P2 


(2.7) 
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| xCo2  = 137  .339  + 0 . 441133  T 

+ .001129  87  T2  + 0.0628105  P 
+ 0 . 562974  P2 


(2.8) 


\lQi  = 190.395  + 0.650043  T 

- . 000897  542  T2  + 0 . 13  0897  P 
+ . 00132418  P2 


(2.9) 


in  which  /x  is  in  ngr/ (c/n  • s)  , T is  in  °C  and  P is  in 
atmospheres  (Jones,  1984) . 

Chapman  and  Enskog  developed  the  following  equation  to 
predict  the  viscosity  of  multi-component  system 


Here  n is  the  total  number  of  species,  Xj  and  Xj  are  the  mole 
fractions  of  the  species  i and  j . yu;  and  fj.i  are  the  dynamic 
viscosities  at  the  pressure  and  temperature  of  the  system, 
and  M|  and  Mj  are  the  molecular  weights  of  the  species.  This 
equation  has  been  shown  to  produce  results  within  2%  of  /xmix 
values  (Bird,  et  al.  1960) . 


(2.10) 


in  which 


(2.11) 


11 


The  kinematic  viscosity  can  be  calculated  from  v = -£, 

P 

where  for  conditions  sufficiently  far  removed  from  the 
critical  point,  and  from  phase  transitions,  the  equation  of 

p 

state  p = — may  be  used. 
r RT  1 


Table  2.3  Predicted  and  Reference  Viscosity 


Gas 

P eq 2.5 

10  ~4g 
cm  • sec 

P 

K RT 
kg 
m 3 

V eq2 . 5 

cm2 

sec 

Pref 
10  -4g 
cm  • sec 

Pref 

kg 

m 3 

vref 

cm2 

sec 

Air 

1.837 

1.172 

0.1567 

1.843 

1.177 

0.1566 

C02 

1.524 

1.789 

0.0852 

1.520 

1.788 

0.085 

02 

2 . 057 

1.301 

0.1581 

2.059 

1.299 

0.1585 

n2 

1.767 

1.139 

0.1551 

1.767 

1.138 

0.1553 

pref  (Edwards  et  at .)  vref  (Johnson,  McCloskey)  p.ref  = Pref  v ref 
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Table  2.4  Diffusion  Collision  Integral  (Skelland  1974) 


kT/e 

kT/e 

kT/e 

0.30 

2 . 662 

1.65 

1.153 

4 . 0 

0.8836 

0.35 

2.476 

1.70 

1.140 

4 . 1 

0.8788 

0.40 

2 .318 

1.75 

1.128 

4 . 2 

0 . 8740 

0.45 

2 . 184 

1.80 

1.116 

4 . 3 

0.8694 

0.50 

2 . 066 

1.85 

1.105 

4.4 

0.8652 

0.55 

1.966 

1.90 

1.094 

4 . 5 

0.8610 

0.60 

1.877 

1.95 

1.084 

4.6 

0.8568 

0.65 

1.798 

2 . 00 

1.075 

4 . 7 

0.8530 

0.70 

1.729 

2 . 1 

1.057 

4.8 

0.8492 

0.75 

1.667 

2 . 2 

1.041 

4.9 

0.8456 

0.80 

1.612 

2 . 3 

1.026 

5.0 

0.8422 

0.85 

1.562 

2.4 

1.012 

6 

0.8124 

0.90 

1.517 

2.5 

0.9996 

7 

0.7896 

0.95 

1.476 

2 . 6 

0.9878 

8 

0.7712 

1.00 

1.439 

2.7 

0.9770 

9 

0.7556 

1.05 

1.406 

2 . 8 

0.9672 

10 

0.7424 

1.10 

1.375 

2.9 

0.9576 

20 

0.6640 

1.15 

1.346 

3 . 0 

0.9490 

30 

0.6232 

1.20 

1.320 

3 . 1 

0.9406 

40 

0.5960 

1.25 

1.296 

3.2 

0.9328 

50 

0.5756 

1.30 

1.273 

3 . 3 

0.9256 

60 

0.5596 

1.35 

1.253 

3.4 

0.9186 

70 

0.5464 

1.40 

1.233 

3.5 

0.9120 

80 

0.5352 

1.45 

1.215 

3 . 6 

0.9058 

90 

0.5256 

1.50 

1.198 

3 . 7 

0.8998 

100 

0.5170 

1.55 

1.182 

3 . 8 

0.8942 

200 

0.4644 

1.60 

1.167 

3.9 

0.8888 

400 

0.4170 
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Table  2.5  Viscosity  Collision  Integral  (Bird,  et  al.  1960) 


kT/e 

kT/e 

kT/e 

0.30 

2.785 

1.65 

1.264 

4 . 0 

0.9700 

0.35 

2 . 628 

1.70 

1.248 

4 . 1 

0.9649 

0.40 

2.492 

1.75 

1.234 

4 . 2 

0.9600 

0.45 

2 . 368 

1.80 

1.221 

4 . 3 

0.9553 

0.50 

2 . 257 

1.85 

1.209 

4.4 

0.9507 

0.55 

2.156 

1.90 

1.197 

4.5 

0.9464 

0.60 

2 . 065 

1.95 

1.186 

4 . 6 

0.9422 

0.65 

1.982 

2 . 00 

1.175 

4.7 

0.9382 

0.70 

1.908 

2 . 1 

1.156 

4.8 

0.9343 

0.75 

1.841 

2 . 2 

1.138 

4 . 9 

0.9305 

0.80 

1.780 

2.3 

1.122 

5.0 

0.9269 

0.85 

1.725 

2.4 

1.107 

6 

0.8963 

0.90 

1.675 

2.5 

1.093 

7 

0.8727 

0.95 

1.629 

2.6 

1.081 

8 

0.8538 

1.00 

1.587 

2 . 7 

1.069 

9 

0.8379 

1.05 

1.549 

2 . 8 

1.058 

10 

0.8242 

1.10 

1.514 

2.9 

1.048 

20 

0.7432 

1.15 

1.482 

3 . 0 

1.039 

30 

0.7005 

1.20 

1.452 

3 . 1 

1.030 

40 

0.6718 

1.25 

1.424 

3.2 

1.022 

50 

0.6504 

1.30 

1.399 

3 . 3 

1.014 

60 

0.6335 

1.35 

1.375 

3.4 

1.007 

70 

0.6194 

1.40 

1.353 

3 . 5 

0.9999 

80 

0.6076 

1.45 

1.333 

3 . 6 

0.9932 

90 

0.5973 

1.50 

1.314 

3.7 

0.9870 

100 

0.5882 

1.55 

1.296 

3.8 

0.9811 

200 

1.60 

1.279 

3 . 9 

0.9755 

400 

CHAPTER  3 


ANALYTICAL  SOLUTION  FOR  ENHANCED 
DISPERSION  IN  OSCILLATING  FLOW 

The  molecular  diffusion  coefficient  and  kinematic 
viscosity  are  two  properties  that  will  arise  in  connection 
with  the  analytical  solution  to  the  enhanced  axial 
dispersion  in  tubes  to  be  developed  next.  Their  values  will 
be  directly  determined  from  the  relation  to  be  devised 
below. 

Velocity  Profile  Solution 

Oscillating  flow  of  a fluid  with  a axial  concentration 
gradient  is  controlled  by  both  the  Navier-Stokes  Equation 
and  Fick's  Law.  The  Navier-Stokes  Equations  control  the 
physical  flow  of  the  fluid  as  a whole,  while  Fick's  Law 
controls  the  diffusion  within  this  fluid. 

The  Navier-Stokes  equation  for  an  incompressible  fluid 
is 

= - -Vp  + v V2u  (3.1) 

Dt  p v ' 
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where 


D 

Dt 
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is  the  material  derivative,  ^ is  the  gradient 


operator,  and  V2  is  the  Laplacian  operator.  In  cylindrical 
coordinates  the  Navier-Stokes  equation  reads 


where  the  velocity  u = erur  + eeue  + exux  and  er , e0 , ex  are 

unit  base  vectors  in  the  radial,  angular,  and  axial 
direction,  respectively.  For  pure  oscillating  flow  confined 
to  a tube  of  infinite  length  and  circular  cross  section,  the 
fluid  velocity  is  independent  of  angular  direction  and  axial 
location.  These  assumptions  reduce  the  Navier-Stokes 
equation  to 


Here,  u = ux  is  the  velocity  in  the  x direction  along  the 


(3.3) 


length  of  the  pipe.  If  it  is  assumed  that  the  pressure 
gradient  varies  periodically  in  time  as 
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1 dP 
p dx 


Kei,J>t 


(3.4) 


where  K = 


1 dP 
p dX 


is  a constant  and  a measure  of  the  magnitude 


of  the  pressure  gradient,  and  the  velocity  is  taken  to  vary 
as 


u (t)  , t)  = $t[£/0  f (rj)  ei45t]  , 


(3.5) 


z , . . 

where  ri  = — is  a non-dimensional  radius.  With  these 
a 


assumptions,  equation  3.3  becomes 

* ^ f = K + v2/ 
v C70v 


(3.6) 


or 

d2f  M 1 df  _ iua2  f __  Ka2  . 

dr\2  t!  5t!  v tf0v  • * } 

An  analytical  solution  can  be  found  for  this  equation  using 
the  appropriate  boundary  conditions  of  f{  0)  = finite  and 
f(  1)  = 0 (no  slip  condition)  (Schlichting  1960) . 
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The  solution  of  equation  (3.7)  is 


f (rj ) = A J0 


"N 


iua‘ 


+ B Yn 


~ i 0)  a 


K 


U0U 


(3.8) 


Introducing  the  non-dimensional  Womersley  Number 


and  applying  the  above  boundary  conditions,  one  finds 


f(ri) 


- i K 
Uq  w 


/ 

1 - 
V 


Jo^Tari)' 

J0(\Fia)  , 


(3.9) 


or  since 


u(ri,t)  = 3t[U0f{ ti)eiwt],  that 


at 

- 

' ± _ J0(v^aT l)  "l 

0) 

1 ) 

(3 . 10) 


Generally,  Jn(/-Ix) can  also  be  written  as  the  Kelvin 
functions  bern{x) +i  bein(x)  where  bern{x)  is  the  real  portion 
of  the  function  and  bein(x)  is  the  imaginary  part.  One  can 
then  rewrite  equation  (3.9)  and  (3.10)  as 


j -iK  _ ber0(gx])  + ibei0  (ar\) 


Uq  w 


bern  (a)  + i beia  (a) 


(3 . 11) 


so  that 
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u (ti  , t)  = - i — 

(0 


her0(ccr| ) + i bei0(ur]  ) 
ber0  (a  ) + i bei0  ( a ) 


(3.12) 


Figures  3.1  and  3.2  show  log  plots  of  and 

JjtyRx)  which  have  a general  shape  of  a spiral.  The 

functions  deviate  from  the  spiral  shape  only  when  the 
independent  variable  is  below  a value  of  about  x = 5.  This 
prediction  of  the  Kelvin  functions  has  heretofore  not  been 
observed,  but  follows  directly  from  the  asymptotic  form  of 
bern{x)  and  bein{x)  as  given  in  Abramowitz  and  Stegan 

(1972).  For  large  x one  finds 


(3 . 13) 
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Figure  3.1  Plot  of  Log  J0( (-i) 5 x) 

90* 


(Numbered  Points  are  values  of  x) 
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Figure  3.2  Plot  of  Log  J,((-i)  5 x) 


270' 


(Numbered  Points  are  values  of  x) 
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Figure  3.3  has  the  various  Kelvin  functions  plotted  for 
the  independent  variable  from  zero  to  six.  This  is  the 
range  of  variables  where  figures  3.1  and  3.2  deviate  from 
the  spiral  shape. 

Figure  3.3  Plot  of  ber0,  bei0,  ber1(  and  beij 


0 1 2 3 4 5 6 


x 

Figures  3.4  through  3.9  are  velocity  profiles  given  by 
equation  3.12  at  various  specified  phases  of  the  pressure 
gradient.  Note  that  the  boundary  layer  thickness  is  reduced 
at  higher  Womersley  Numbers. 
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Figure  3.4  Velocity 
Profile  at  0 Degree 


Profile  at  60  Degree 


Profiles  at  120  Degree 


Figure  3.5  Velocity 
Profile  at  30  Degree 


Profile  at  90  Degree 


Profile  at  150  Degree 


Velocity  Profiles  at  Specified 
Pressure  Phase  for  Various 
Womersley  Numbers 
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Concentration  Solution 


Using  this  velocity  profile  in  Fick's  Equation  for 
diffusion 


d£ 

dt 


+ 


(U7S) c * °»V2C 


(3 


and  making  the  approximation  for  concentration  of 

C( i\,x,  t)  = y(x  + Eg(r\)ei(jC)  (3 


which  was  introduced  by  Chatwin  (1975),  equation  3.14 
becomes 


i (0 


9 + 


Un 


f = V2gr 


(3 


or 


d2g  + _ldg  _ itoa2  _ a2uo  ^ 
6ti 2 n 3ti  dm  dme 


Note  that 


to  a‘ 


D, 


can  be  re-written  as 


w v _ „2 


M 


v D, 


- a S where 


M 


the  Schmidt  number,  so  that  one  has  the  alternate  form 


&9  + ±M-ia*Sg  = £Elf 

dr| 2 r\  at)  DmE 


.14) 


. 15) 


. 16) 


.17) 


S is 


(3.18) 
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This  equation  is  subject  to  the  boundary  conditions 
= 0 (no  diffusion  through  the  wall)  and 

finite.  Assuming  a solution  of  the  form 

g(  T|)  = (f(Ti)  + h(r\))F  , (3.19) 


dgr(l) 

aii 

g{o)  = 


equation  3.18  becomes 


d{f  + h)  + 1 d{f  + h) 

dr\2  r)  9ti 


ia2S{f  + h) 


*2Uq  f 

DmEF 


(3.20) 


Using  equation  (3.7),  the  partial's  of  f can  be  eliminated 
resulting  in 


d2h  + _1  dh 

dr]2  t)  dr) 


i a2  Sh 


a2U, 

DmEF 


+ i a2  (5-1) 


f + 


Fa2 

C70v 


(3.21) 


Setting  F = i the  f term  drops  out  leaving  us 

E co  ( 5 - 1 ) 


with  the  equation 


d2h  + _ldh 

at)  rj  dt] 


ia2Sh 


Ka2 
Uq  v 


(3.22) 
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The  solution  to  this  equation  is 


h{i\) 


= GJ0  {^-iS a t])  + HY0  (yZ-iSa  tj)  +i 


K 

£70  co  S 


(3.23) 


To  satisfy  g{  0)  = finite  the  parameter  H must  equal  zero. 
Also  from  the  Bessel  function  derivative  identity 
dJQ  ( x ) 


dx 


- Jx  (x)  , we  have 


dJ0  (y/-iSa t]) 

a^T 


- y/-iS  a Jx  (\f-TS  at) ) 


(3.24) 


The  derivative  for  g{r\)  is 


dg{r\) 

di] 


= 0 = 


T)  =1 


u0s 


E co  ( S - 1] 


- G\J-i  S a Jx  (^-i  S a) 


J0  a) 


(3.25) 


Solvinq  for  G,  one  finds 


KJ^y/^ia) 


UQ  co  v/5  Jq  (v^-l  a ) Jx  {^-iS  a ) 


(3.26) 
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Therefore 


^(»l) 


K 


\[3  a)  J0W~i  Su  r|)  _ S J0(/=i  a x\) 


\ 


(3.27) 


£u2  (S  - 1) 


J0  (vFI  a)  J-!  W-iS a)  J0  a) 


Determining  the  Enhanced  Diffusion  Coefficient 


The  enhanced  coefficient  of  diffusion  for  oscillatory 
flow  can  be  determined  by  setting  the  effective  diffusion 
equal  to  the  time  averaged  net  concentration  flux  through 
the  cross-sectional  area,  of  the  tube.  Mathematically  one 
has 


(3.28) 


In  this  integral,  u = UQ  {uR  + i ux)  , where  uR 


feiat  + fe~iut 
2 


and  Uj 


f ei<Jt  - fe‘iwt 


, and  C = CR  + i Cj , where 


2 


2 
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The  bar  denotes  the  complex  conjugate  of  the  function  and  R 
and  I refer  to  the  Real  and  Imaginary  components 
respectively.  Equation  (3.28)  may  thus  be  re-written  as 


n z dC  _ TT  r I f eiut  + f e-iut  ^ . feiuc  - fe-iwt\ 

D°ffA  " U°  J \ 2 1 2 ) 

|z  + ggeiut  + gge'iMt:|  + i Eg ei<Jt  - Eg 


(3.29) 


or 


DeffA^  = U0  y (1  felat*  (1  ~i} 


+ -LL.+  ..4 ) e g eiu  ^ ^ il. 


E’gre 


-ito  tj 


. (3.30) 


dA 


When  integrated  over  the  one  period  of  oscillation  to  obtain 
a time  averaged  value,  the  remaining  terms  will  be  time 
independent.  The  resulting  time  averaged  integral  reads 

Deff  A = E U°  Y J (f  (i))  ^(ti)  + f[i\)g{i\))dA  (3.31) 


The  application  of  Green's  Theorem  will  transform  this 
last  integral  into  a form  that  is  more  easily  solved.  When 
equation  (3.6)  is  multiplied  by  the  conjugate  of  g and  the 
conjugate  of  equation  (3.16)  is  multiplied  by  f,  one 
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finds 


Eg  = _2L  g + gV2f 
v a U0v 


(3.32) 


and 


- _L!f£  f g+  f = jfV2^  . 

\ DE 


(3.33) 


When  (3.33)  is  subtracted  from  (3.32),  one  has 


) fg-  ff  = g + gV2f  - fV2g  . (3.34) 

Dm  v Dwi?  v t70  y 


It  is  desirable  to  eliminate  g from  this  equation.  This 
can  be  done  by  inserting  into  equation  (3.16)  the  equivalent 
of  f from  equation  (3.6),  taking  the  conjugate  of  the 
result.  That  is, 


KU0 

co2  E 


v £7n  i D, 


CO 


2 E 


V2f  + *S/2g 


co 


(3.35) 
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This  yields 


gV2f  - f V2g  + .i5.ffgvay+  — ^-V2f 


w v C7n 


co2£ 


(3 


LG)  x f *T2 


, 1 CO  ICO  > — 

( + _±lt  ) f g - ^ ^ 

v A/  dme  \ <ss2  E 


and  the  conjugate  of  the  above  is 


l KD, 


gV2f-fV2g-  J"'~MV2g+  — *-V2f  = 


CO  V Ur 


CO2  E 


(3 


- (A^  + 4^)f  g-ELff- 


K2 


D 


M 


dme 


\ (j32  E 


Now  (3.37)  is  subtracted  from 


2 n 

g^f  - g^f  - fV 2g  + f V2^  + 

^ CO  V CL 


1 CO  + 1 CO 


(3.36)  which  gives 

('V2g+'V2g)  + — E—  (V2/ " V2/) 
CO2  E 

(fg+  fg) 


.36) 


. 37) 


(3.38) 
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Equation  3.38  is  inserted  into  equation  3.31  which  yields 


DeffA 


u0 dm se 


dc 


dz  i 0)  (5  + 1)  2 


/( 


g^f  - gV*f  - f^g  + f V2^  + KDm  (V2^  + V2^)  + — (V2^  - V2/)  |dA 


w v f/n 


<o2  E 


(3 . 39) 


Utilizing  Green's  Second  Theorem  gives 


dc 


U0DmSE  _ 


eff  dz  l co  (5+1)4 


<f(  gVf-gVf-fVg  + fVg+  (Vg  + Vg)  + — (Vf  - V/)  Ids 

J l <j)  v Un  a2  E 


(3.40) 


And  using  the  boundary  condition 


at  the  tube  wall 


yields 


■'eff 


U°D»SE  Yff/j- 

i co  (5+  1)  4 J |\a 


— 1 
co2  E) 


Vf - 


(3.41) 
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If  two  new  functions  are  defined  as 


G{  tj)  =l-^g[r\) 


and 


F(i\)  = l^-f(  T|) 
K. 


(Watson  1983)  the  integral  is  reduced  in 


complexity  to 


D 


eff 


Dm.SK  Y <f(GV  F + GV  F)  dS 
4 (5  + 1)  co4  • 


(3.42) 


where  dS  = T]  d0 . This  equation  is  solved  simply  as 


ac  . _E»££!UL(g(i)  W ii 

dx  2 (5  + 1)  co4  \ 3ti  3ri 


(3.43) 


where 


and 


dF(l)  = _ ■ m U0  df(  1) 
5t|  K dt| 


It  is 


interesting  to  note  that  since  this  is  the  sum  of  a function 
and  its  conjugate,  the  result  is  the  real  portion  only. 


1 

f -VSJ1(Pia)J0(PTSa) 

(5-1) 

{ Jq  ip*  a ) J\  W~i  S oc ) 

(3.44) 


and 


dF{l)  = _ P « a) 

3tl  ~ a) 


(3.45) 
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Multiplying  equation  3.44  and  equation  3.45  produces 


G(i)  = 

dr] 


1 

- \/S  Jx  (\/-ri  a ) J0  (v/-iSa) 

' sfi  a J1(v/772'  a)N 

(1  - S) 

{ J0W~i  a)  Ji  W~i  Sa) 

k Jo(v^ «)  , 

Taking  twice  the  real  portion  of  equation  (3.46)  will 
provide  the  sum  of  the  function  and  its  conjugate  of 
equation  (3.43),  namely 
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Tidal  Displacement  Solution 

The  average  tidal  displacement  can  be  calculated  by 
taking  the  magnitude  of  the  integral  of  one-half  of  the 
period  times  the  area  of  the  pipe  and  equating  it  to  the 
tidal  displacement  times  the  area  of  the  pipe.  We  have 
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and 
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Due  to  the  fact  that  Jxni7n_1  (x)  dx  = xnJn(x)  (Abramowitz  and 


Stegan) , upon  integrating  we  find  (Kurzweg  1988) 


2 K 
0)2 


1- 


2 i 

v^ia  I d0  a) 


Ax  = 


(3.51) 


34 


or 
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A plot  of  equation  (3.52)  relating  axial  pressure 
gradient  K to  tidal  displacement  is  shown  in  Figure  3.10  and 
3.11.  Figure  3.12  is  a phase  relationship  between  the 
pressure  gradient  and  velocity  represented  by 
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Figure  3.10  Pressure-Displacement  Relationship  (Smaller  a) 
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Figure  3.11  Pressure-Displacement  Relationship  (Larger  a) 
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Phase  of  Pressure  - Velocity 
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Enhanced  Dispersion  Coefficient  in  Final  Form 


Defining  the  non-dimensional  dispersion  coefficient 
X = Deff/ co Ax2,  we  have 


X [a , S)  = 
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Now,  letting 


F(a)  = F*(a)  + i Ft [a) 


^ Berg  (a)  +iBeio(u)  _ . Joi/^a)  _ _ . /7iJ1(vci«) 
2 Ber0  ( a ) + iBei0(oc)  1 ^(/^a)  2 J0  (/^T  aP 


equation  (3.54)  may  be  re-written  in  it's  final  desired  form 
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This  version  for  effective  or  enhanced  dispersion 
coefficient  in  oscillatory  flow  was  first  obtained  by 
Kurzweg  and  Jaeger  (1987) . 
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Graphical  Presentation 


Figure  3.13  is  a contour  plots  of  A = Deff/ to  Ax  as  a 
function  of  Womersley  Number  a and  Schmidt  Number  5.  Two 


curves  = 0 and  44  = 0 are  plotted  as  well.  4^-  = 0 

0a  05  ^ 0a 


specifies  the  maximum  value  of  X for  a given  Schmidt  number. 

dX 

-5-  = 0 specifies  the  maximum  X for  a given  Womersley 
dS 


Number.  These  curves  are  useful  for  studies  involving  the 
separation  of  species,  showing,  for  example,  that  no 

diffusional  separation  will  be  possible  when  3-  = 0. 


Schmidt 
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Figure  3.13  Contour  Plot  of  X = Deff/co  Ax2 


CHAPTER  4 


EXPERIMENTAL  PROCEDURE 

In  the  previous  Chapter  we  have  derived  the  theoretical 
enhanced  dispersion  coefficient  in  oscillatory  flow.  This 
result  will  now  be  tested  experimentally. 

Enhanced  Diffusion  Process 


Enhanced  diffusion  is  a process  by  which  a species  in  a 
fluid  oscillating  in  a tube  can  experience  a vastly 
amplified  rate  of  axial  dispersion  as  compared  to  that 
provided  by  normal  molecular  diffusion.  This  enhancement  is 
promoted  by  an  interaction  between  the  boundary  layer  and 
the  core  as  the  fluid  oscillates.  The  species  under 
consideration  is  transferred  first  by  diffusing  radially 
into  the  boundary  layer  and  later  in  the  cycle  diffusing, 
again  radially,  back  into  the  oscillating  core  of  the  fluid, 
being  transported  axially  a large  distance  from  the 
originating  location. 

The  enhancement  is  caused  through  an  effective  increase 
in  diffusing  surface  area  through  the  stretching  action  due 
to  shear  in  the  fluid,  the  boundary  layer  thickness,  time 
allowed  for  diffusion  to  occur  each  cycle,  and  to  a very 
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limited  degree  an  increase  in  the  concentration  gradient. 
This  increase  in  concentration  gradient  can  be  visualized  by 
considering,  hypothetically,  a tube  with  a linear 
concentration  variation  of  one-hundred  percent  per  one- 
hundred  centimeter  in  a tube  of  two  centimeter  diameter.  If 
the  tidal  displacement  caused  the  core,  which  will  remain  at 
zero  percent,  to  approach  the  area  where  the  fluid  at  the 
tube  wall  will  remain  at  one-hundred  percent  due  to  the  no- 
slip requirement,  the  overall  radial  concentration  gradient 
will  now  be  one-hundred  percent  per  centimeter. 

The  flow  in  an  oscillating  fluid  is  influenced  by  the 
pressure  gradient,  inertial  forces,  and  viscous  forces.  The 
boundary  layer  thickness  is  determined  by  the  frequency  of 
oscillation  and  the  kinematic  viscosity  of  the  fluid.  As 
the  frequency  is  increased,  the  boundary  layer  thickness  is 
decreased.  Likewise,  the  boundary  layer  thickness  is 
decreased  as  the  kinematic  viscosity  is  decreased.  This 
thickness  influences  the  total  amount  of  a species  that  it 
can  contain.  As  the  thickness  of  the  boundary  layer 
decreases,  the  capacity  to  hold  a species  is  reduced  due  to 
the  reduced  mass  of  the  boundary  layer.  The  amount  that  the 
boundary  layer  will  contain  is  also  limited  by  the  amount 
that  a species  that  can  diffuse  into  the  boundary  layer. 
There  is  a diffusion  layer  thickness  that  will  contain  the 
species  that  has  diffused  into  the  layer.  If  the  boundary 
layer  thickness  is  much  greater  than  the  diffusion 
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thickness,  little  mass  of  the  species  will  be  contained  in 
the  boundary  layer.  Likewise,  if  the  boundary  layer  is  much 
less  than  the  diffusion  thickness,  little  mass  will  be 
transferred,  due  to  the  restriction  of  the  boundary  layer 
thickness.  Therefore,  there  is  a maximum  mass  transfer  when 
the  boundary  layer  thickness  approximately  matches  the 
diffusion  thickness.  The  amount  of  mass  of  a species 
transferred  to  the  boundary  layer  is  also  influenced  through 
the  time  allocated  for  transfer  during  a cycle.  When  the 
oscillating  frequency  is  small,  a larger  time  is  allowed  for 
a species  to  transfer  to  the  boundary  layer.  This  is  one 
factor  which  determines  the  diffusion  thickness.  As  the 
frequency  is  increased,  the  boundary  layer  thickness 
decreases,  as  well  as  the  diffusion  thickness  decreasing  due 
to  the  reduced  time  allowed  for  diffusion.  The  above 
peaking  of  the  lambda  curve  (figure  3.13)  is  thus  determined 
by  the  matching  of  boundary  layer  thickness  to  the  diffusion 
thickness  and  leads  to  tuning  effect  first  discussed  by 
Kurzweg  in  1986. 

The  pressure-velocity  phase  relationship  (figure  3.12) 
is  due  to  the  relative  change  in  importance  of  inertia  and 
viscous  forces.  The  ratio  of  inertial  to  viscous  forces  in 
oscillating  pipe  flow  goes  as  the  square  of  the  Womersley 
number.  The  viscous  forces  are  greater  than  inertial  forces 
when  the  frequency  is  low.  As  the  frequency  is  increased, 
the  inertial  forces  increase  faster  than  viscous  forces, 
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eventually  overtaking  them.  These  two  forces  are  inherently 
90°  out  of  phase,  therefore,  the  phase  change  is  simply  a 
result  of  the  changing  importance  of  the  two  forces. 

The  Pressure-Displacement  curve  (figure  3.10  and  3.12) 
is  also  a result  of  the  phase  change  between  the  viscous 
forces  and  inertial  forces.  This  curve  is  experimentally 
determined  by  measuring  the  magnitude  of  the  pressure 
gradient,  and  the  tidal  displacement  Ax  of  the  fluid. 

Along  with  the  fluid  density,  frequency  of  oscillation  and 
tube  radius,  a curve  can  be  generated.  This  curve  is 
independent  of  the  resonance  of  the  system,  therefore,  it 
can  be  used  to  accurately  infer  the  tidal  displacement. 


Experimental  Procedure 


The  experimental  procedure  was  to  measure  the  enhanced 
diffusional  flux,  and  calculate  the  associated  enhanced 
axial  dispersion  coefficient.  The  diffusional  flux  was 
measured  in  two  separate  setups  with  frequency,  tube  radius, 
and  tidal  displacement  varied  over  a wide  range  of  values. 
The  imposed  constraints  were  to  limit  the  tidal  displacement 
to  a maximum  of  one-third  of  the  length  of  the  tubes,  and 
conditions  set  to  remain  in  the  laminar  flow  regime.  To 
remain  in  the  laminar  regime  a non-dimensional  term 

— < 80  0 had  to  be  met  at  all  times.  The  turbulent 
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range  was  determined  by  Kurzweg,  Lindgren  and  Lothrop  in 
1989. 

The  diffusion  flux  rate  was  determined  by  measuring  the 
concentration  in  the  reservoirs  with  the  associated  flow 
rates.  The  measured  diffusion  coefficient  can  be  obtained 
from  the  diffusional  flux  rate,  cross-sectional  area,  length 
of  tube,  and  the  tidal  displacement  either  measured 
directly,  or  from  the  measured  magnitude  of  pressure 
gradient. 

The  experiment  was  performed  on  two  different 
experimental  setups  shown  in  figures  4.1  through  4.3.  The 
first  used  nitrogen  from  the  atmosphere  diffusing  into 
oxygen.  This  apparatus  was  used  to  acquire  initial 
operating  experience  with  the  equipment  and  for  preliminary 
trial  runs.  The  second,  more  sophisticated  apparatus  used 
by  Dr.  Jaeger  for  gas  separation  studies,  was  used  here  to 
obtain  data  on  oxygen  diffusing  into  carbon  dioxide. 

The  first  setup  used  a measure  of  pressure  gradient 
instead  of  tidal  displacement.  This  apparatus  consisted  of 
three  different  sized  PVC  tubes  with  two  taps.  These  taps 
were  located  near  the  center  of  the  tube  in  an  attempt  to 
reduce  the  influence  of  'end  effects'  that  consists  of  three 
major  components.  These  are  entrance  turbulence,  non  fully 
formed  velocity  profiles  and  possible  concentration  gradient 
perturbations.  The  expected  side  effects  of  center  taps 
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along  the  length  are  reduced  pressure  differences  and 
concentration  differentials. 

Figure  4.1  Schematic  of  Setup  1 
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The  second  experimental  apparatus  consisted  of  a system 
of  76  glass  capillaries,  each  being  0.38  centimeters  in 
diameter  and  120  cm  in  length.  The  tidal  displacement  was 
measured  by  means  of  a spirograph  which  had  an  electrical 
output  proportional  to  the  volumetric  displacement.  The 
measurement  of  tidal  displacement  was  made  independent  of 


Q3 
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Figure  4.2  Schematic  of  Setup  2 


the  actual  diffusion  test  runs,  as  it  was  felt  that  the 
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additional  load  on  the  system  would  affect  the  sinusoidal 
wave  form  and  cause  possible  compressibility  effects.  It 

Figure  4.3  Photograph  of  setup  2 


was  not  possible  to  obtain  a true  pressure  gradient  with  the 
second  setup  due  to  the  tube  bundle  construction,  however, 
pressure  was  measured  in  the  reservoirs  and  end  effects  were 
assumed  to  be  minimal.  Since  the  system  had  to  be  run  as  a 
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closed  system  with  oxygen  diffusing  into  carbon  dioxide,  a 
more  complex  experimental  apparatus  was  involved. 
Measurements  of  the  inflow  and  outflow  for  each  end 
reservoirs  had  to  be  measured  and  controlled. 

EXPERIMENTAL  ASSUMPTIONS 


Diffusion  Coefficient  is  Independent  of  Concentration 

The  Diffusion  Coefficient  does  not  change  appreciably 
with  concentration  changes,  but  a minor  change  does  occur 
(Bird  1960)  (Hirschf elder , et  al.  1949).  Since  this  is  a 
minor  change,  an  assumption  of  constant  molecular  diffusion 
coefficient  will  not  have  a great  influence  in  the  final 
results . 

Incompressible  gases 

The  assumption  that  one  is  dealing  with  an 
incompressible  gas  was  an  assumption  that  was  not  quite  met 
in  this  experiment,  but  closely  approximated.  The  worst 
departure  from  the  incompressible  gas  assumption  occurs  for 
the  high  tidal  displacement  and  high  frequency  operation,  as 
these  are  the  conditions  of  greatest  pressure  gradient 
requirements.  The  great  majority  of  the  experimental  runs 
had  a pressure  that  fluctuated  at  most  by  0.04  atmospheres. 
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Such  conditions  can  be  assumed  to  meet  the  incompressible 
gas  condition. 

The  effects  of  compressibility  include  velocity  changes 
due  to  the  compression  of  the  gas,  the  change  in  the 
molecular  diffusion  coefficient  due  to  pressure  changes, 
which  will  affect  the  Schmidt  number. 

Linear  Concentration  Gradient 


Although  the  concentration  gradient  will  vary  in  tubes 
where  a counter flow  exists,  the  concentration  gradient  can 
be  considered  as  linear  for  a short  length  from  any  given 
plane  normal  to  the  tube  axis.  The  actual  concentration 
gradient  in  the  first  setup  is  an  exponential  curve  along 
the  length  of  the  tube.  In  the  second  setup,  the 
concentration  gradient  is  in  fact  linear.  This  assumption 
is  best  approximated  when  the  tidal  displacement  is  small. 

Pressure  Gradient  is  a Function  of  a Pure  Sine  Wave  in  Time 


A large  problem  in  this  experiment  was  in  generating 
pure  sine  waves.  A linear  motor  was  used  to  generate  a pure 
driving  sine  wave,  shown  in  figure  4.4,  but  there  were  two 
separate  causes  of  harmonics  in  this  experiment.  One  is  the 
acoustic  freguency  of  the  tube  and  the  other  is  the  tube  and 
isolator  coupling. 
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Figure  4.4  Photograph  of  Linear  Motor  Used  to  Produce 
Sinusoidal  Axial  Displacements  of  Fluid  within  Conduit. 


The  acoustic  frequency  of  the  tube  is  caused  by  a one- 
quarter  wavelength,  or  an  odd  multiple  of  a quarter 
wavelength,  pressure  wave  standing  in  the  tube.  This  is 
controlled  largely  by  the  physical  length  of  the  tube. 
Lengthening  of  the  tube  causes  a lowering  of  the  acoustic 
frequency.  Mathematically,  the  acoustic  frequency  , i.e. 
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Figure  4.5  Photograph  of  Pressure  Variation  (Top) 
and  Tidal  Volume  (Bottom) 


0.5  Second 


the  quarter  wave,  is  w = — - , where  c is  the  acoustic 

4 L 

velocity  and  L is  the  length  of  the  tube. 

Another  cause  of  harmonics  is  from  a coupling  of  the 
tube  and  the  end  chambers.  This  frequency  is  caused  from 
the  compressibility  of  the  gas  in  the  end  chamber  acting  on 
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the  mass  of  the  gas  oscillating  in  the  tube.  In  the  non- 
viscous  case,  this  freguency  is  directly  proportional  to  the 
diameter  of  the  tube  and  inversely  proportional  to  the 
square  root  of  both  the  length  of  the  tube  and  the  volume  of 
the  isolation  chamber.  The  non-viscous  natural  frequency  is 

(o  = a C where  a is  the  radius  of  the  tube,  c is  the  speed 

v/1 7V 

of  sound,  L is  the  length  of  the  tube  and  V is  the  volume  of 
the  isolation  chamber. 

The  acoustic  frequency  provided  high  frequency 
harmonics  generally  of  the  tenth  order  that  caused  the 
driving  frequency  to  have  imposed  upon  it  this  tenth  order 
harmonic.  (See  Figure  4.5)  The  tube-isolation  chamber 
frequency  was  of  the  same  order  of  frequency  as  the  driving 
frequency,  the  main  problem  being  a resonance  that  causes 
the  tidal  displacement  to  be  much  larger  than  it  otherwise 
would  have  been.  Since  the  pressure  gradient  is  related  to 
the  tidal  displacement,  no  adverse  effect  was  observed  by 
this  when  using  the  pressure-tidal  displacement  measurement, 
provided  the  waveform  remained  as  a pure  sine  wave. 

The  acoustic  frequency  was  the  larger  of  the  two 
problems  when  the  driving  frequency  was  high.  The  solution 
was  to  use  a short  length  of  tube,  forcing  the  acoustic 
frequency  to  become  very  high.  The  effects  eventually 
became  negligible  when  a short  tube  was  used.  At  low 
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frequencies,  the  isolator-tube  coupling  was  the  larger 
problem.  There  was  no  solution  to  this  problem  for  the 
equipment  on  hand  except  to  keep  the  operating  frequency 
above  the  point  where  the  harmonics  became  dominant. 

Cross-Sectional  Velocity  Profile  is  Fully  Developed 

Taking  pressure  gradient  measurements  near  the  center 
of  the  tube  length  insures  that  the  velocity  profile  is  as 
fully  developed  as  can  be  obtained  in  this  experiment.  This 
is  one  of  the  'end  effects'  problems  mentioned  earlier. 
Indications  from  earlier  studies  (Kurzweg,  Lindgren  and 
Lathrop  1989)  indicate  end  effects  in  tube  length  only 
extend  about  Ax  distance  along  the  tube  from  their  ends. 

Transverse  Diffusion  Effects  Overwhelms  Axial  Diffusion 


Transverse  diffusion  is  assumed  in  the  theoretical 
development  to  be  the  only  component  of  molecular  diffusion 
involved  in  producing  of  the  enhanced  diffusion.  Since  the 
enhanced  diffusion  is  several  orders  of  magnitude  larger 
than  molecular  diffusion,  this  is  generally  a safe 
assumption.  This  will  not  be  true  in  exceptionally  low 
Womersley  Number  or  extremely  high  Womersley  Number 
conditions,  where  the  lambda  is  very  low  as  is  suggested  by 
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Figure  3.13.  In  these  cases,  the  enhanced  diffusion  may  be 
comparable  to  molecular  diffusion. 

Flow  is  Completely  Laminar 

It  is  assumed  that  the  flow  is  completely  laminar 
everywhere.  The  effects  of  turbulence  are  not  completely 
understood  in  enhanced  dispersion.  Since  self-maintaining 
turbulence  in  oscillating  flow  generally  develops  slowly, 
the  effects  of  some  limited  end  effect  turbulence  are 
possibly  negligible,  but  must  not  be  forgotten  (Kurzweg, 
Lindgren  and  Lothrop  1989) . 

Steady  State  Conditions  have  been  Achieved 

The  measurements  made  in  this  experiment  are  only  valid 
for  steady  state  laminar  conditions.  The  various 
measurements  were  observed  to  ensure  that  steady  state 
conditions  were  attained.  Small  deviations  from  steady 
state  conditions  possibly  will  not  affect  the  outcome 
appreciably.  This  may  not  be  the  case  if  obvious  non-steady 
state  conditions  exist.  In  the  first  setup,  steady  state 
conditions  sometimes  took  as  long  as  six  hours  to  attain, 
and  therefore,  a few  readings  may  have  been  made  before 
steady  state  conditions  had  actually  occurred. 
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Tube  is  Completely  Straight 

The  tube  is  kept  as  straight  as  possible  to  reduce  the 
effects  of  accelerations  cause  by  the  deviations.  These 
accelerations  will  affect  the  flow  characteristics  of  the 
fluid  that  will  invalidate  the  theory  developed  for  this 
experiment.  Previous  investigators  have  found  that  tube 
curvature  increases  the  rate  of  diffusion  (Sharp  et  al. 
1991) . 

All  Measurements  are  Accurate 


The  calibration  of  the  various  pieces  of  equipment  used 
in  this  experiment  were  calibrated  as  follows.  The  flow 
rate  of  the  feed  gas  was  measured  with  a flow  meter  that  was 
calibrated  by  measuring  the  volume  delivered  over  a 
specified  time.  The  flow  was  calibrated  at  steady 
atmospheric  conditions.  The  frequency  was  checked  at  low 
frequency  and  assumed  to  remain  consistent  as  the  frequency 
was  increased.  The  tube  radius  was  measured  using  a caliper 
and  averaging  the  measurement  over  various  angular 
measurements.  A water  column  was  used  to  calibrate  the 
pressure  transducer.  Although  this  is  a static  measurement, 
it  is  assumed  that  the  calibration  will  remain  valid  for 
transient  measurements.  The  concentration  measurement  could 
be  easily  calibrated  at  zero  and  atmospheric  conditions  only 
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in  the  first  setup,  at  zero  and  one  hundred  percent  for  the 
second  setup.  The  meters  were  assumed  to  be  linear  in  their 
output,  as  no  easy  way  to  calibrate  for  intermediate 
concentrations  was  found. 

No  Leaks  To  or  From  Tube  System 

One  large  problem  early  in  the  experiments,  was  a 
diaphragm  used  in  the  isolation  end  chamber.  This  diaphragm 
was  necessary  to  prevent  atmospheric  gases  from  leaking  in 
or  experimental  gases  from  leaking  out  of  the  reservoir  area 
of  the  isolator  through  the  pump  piston.  Originally, 
Reynold's  Wrap  was  used  as  the  diaphragm  material,  but  was 
found  to  allow  the  gases  to  diffuse  directly  through  this 
diaphragm  material,  invalidating  the  data  collected.  The 
next  diaphragm  used  was  some  thin  rubber  material.  It  was 
found  to  be  too  flexible,  allowing  the  inlet  to  become 
sealed  with  the  diaphragm  material  during  part  of  the  cycle. 
A composite  diaphragm  was  used  for  the  remainder  of  the 
experiment.  This  diaphragm  was  constructed  of  three  layers, 
the  two  outer  layers  were  Reynold's  Wrap  with  the  inner 
layer  being  thin  rubber  material.  This  composite  diaphragm 
met  the  requirements  of  being  impenetrable  to  gas  diffusion 
through  the  material  plus  had  an  appropriate  strength. 


57 


Tidal  Displacement  Remains  Small 

The  tidal  displacement  must  remain  small  compared  to 
the  length  of  the  tube  to  insure  that  flushing  of  the  tube 
is  kept  at  a minimum.  If  the  tidal  displacement  were  to  get 
excessively  large  exceeding  the  length  of  the  tube,  there 
would  be,  in  addition  to  enhanced  diffusion,  considerable 
axial  convection.  To  insure  that  some  unflushed  gases 
remain  in  the  tube,  a limit  of  one-third  the  tube  length  for 
tidal  displacement  was  arbitrarily  chosen  as  the  upper  limit 
used. 


Deriving  the  Non-dimensional  Expressions 
Used  in  the  Experimental  Measurements 
of  the  Two  Setups 

The  first  setup  had  a net  flow  into  the  tube  which 
balanced  the  upstream  diffusion.  This  flow  was  caused  by 
the  injection  of  pure  oxygen  into  the  reservoir.  This  causes 
the  axial  concentration  to  deviate  from  the  usual  linear 
concentration  variation.  A method  has  to  be  developed  to 
get  a numerical  value  for  the  enhanced  diffusion 
coefficient . 

The  diffusion  in  the  first  experiment  is  balanced  by  an 
equal  and  opposite  convective  flow.  Thus  at  mean  velocity 
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U0  one  has 


U0  C 


'eff 


a c 

dz 


(4.1) 


The  general  solution  to  this  equation  is 

U0z 

C = A e °9tt  . 


(4.2) 


To  solve  for  the  Effective  Diffusion  Deff/  the  log  of  the 
above  equation  is  taken.  One  has 

ln(C)  = In  {A)  + U~—  . (4.3) 

Deff 


Boundary  conditions  of  C(0)  = C0  and  C(L)  = Cx  are  applied 
resulting  in 


D. 


eff 


(4.4) 


The  velocity  is  determined  by  the  volumetric  flow  rate  of 
gas  into  the  isolation  chamber 


'eff 


QL 


A In 


/Cl 


(4.5) 


59 


where  Q is  the  volumetric  flow  rate  of  the  gas  into  the 
isolation  chamber,  L is  the  length  between  sample  ports,  A 
is  the  cross-sectional  area  of  the  tube  and  C0 , Cx  are  the 

concentration  measurements  at  the  two  sample  ports.  (See 
Figure  4.1) 

The  Effective  Diffusion  Coefficient  equation  derived  in 
Chapter  3 is  difficult  to  use  with  the  first  experimental 
setup  due  to  the  fact  that  we  are  measuring  the  pressure 
gradient  instead  of  measuring  the  tidal  displacement.  Since 
the  measured  quantity  is  the  pressure  gradient,  it  is  easier 
to  recast  the  equation  into  components  using  the  pressure 
gradient  measurement.  Taking  the  original  function  for  Deff 
given  by  equation  3.47,  namely 


D 


eff 


1 dP 

2 

p dx 

a G)3  (1  - S2) 


J0{J=T5a) 
J0  {/-i  a)  Ji  (V"  i S a) 


\ 

+ s 

) 


, a)  j 


(4.6) 
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a new  function  N can  be  defined  as 


N{a,S) 


gf^l/  -/5J1(/rItt)J0(v'-j  Sa) 
|\  J0  (yFIa)  Jx  (y/-i  Sa) 

4 a (1  - S2) 


\ 

+ 5 

j 


K JoisTia)  J 


(4.7) 
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resulting  in 


N(a,S) 


P2  <*>3  Deff 


4 


dp  2 
dx 


(4.8) 


Substituting  eguation  (4.5)  for  Deff  and  recasting  in 
experimentally  derived  terms,  yields 


N(ct,S) 


QL  p2  cj- 


4 A In 


(CA 

dP 

UJ 

dx 

(4.9) 


Some  reasonable  approximations  can  be  made  to  convert  this 
into  measurable  quantities.  The  pressure  gradient  can  be 


approximated  as 


dP 

dz 


P -P 

M) 


where  L is  the  interval  between 


pressure  ports  in  setup  one  and  a = n a 2 is  the  tube  cross- 


Schmidt 
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sectional  area.  Thus  we  arrive  at  the  desired  non- 
dimensional  expression 


N(a,S) 


0 L3  p2  a)3 

4 71  a2  - P0  )2  In 


(4.10) 


A plot  of  N(a,S)  is  shown  in  Figure  4.6.  It  again  shows 

Figure  4.6 
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a tuning  effect  with  maximum  values  of  N occurring  near 
S = 1 and  a = 4 . 

In  the  second  experimental  setup,  no  net  flow  was  the 
desired  condition,  therefore,  the  diffusion  equation  is 
simply 


= 0 (4.11) 

ax2 


or 


4^  = Constant  . 
ox 


(4.12) 


The  enhanced  diffusion  coefficient  is  therefore  simply 


D. 


eff 


Qx  CXL 
(Co-qjA  ' 


(4.13) 


where  Qx  is  the  low  concentration  outflow,  C0  and  C1  are  the 

high  and  low  concentration  respectively  at  the  capillary 
ends,  L is  the  length  of  the  tubes  and  A is  the  total  tube 
cross  sectional  area. 

The  enhanced  dispersion  equation  (3.55)  is  equated  to 
this  to  obtain  an  experimentally  determined  value  for  lambda 
of 

Q&L 

A (Cg-C-J  Ax2  0) 


X (a , a)  = 


(4.14) 


63 


Determining  Density  and  Viscosity  Experimentally 

It  is  interesting  to  note  that  by  performing  the 
pressure-displacement  experiment,  both  the  density  and 
viscosity  coefficient  of  a gas  under  study  can  be 
determined.  The  process  best  fits  the  collected  data  to  the 
pressure-displacement  curve  with  a least  sum-of-squares  on 
both  parameters  density  and  viscosity,  which  will  provide 
the  density  and  kinematic  viscosity  directly. 

Another  method  utilizes  the  inflection  point  of  the 
pressure-displacement  curve.  The  inflection  point  of  this 
curve  is  at  a Womersley  Number  of  1.62  with  the 
corresponding  function  value  of  0.298.  This  method  would 
best  fit  the  data  to  a quartic  function  that  will  provide  an 
inflection  point.  If  sufficient  data  can  be  collected,  the 
inflection  point  can  be  determined,  assuming  that  the 
pressure  gradient,  frequency,  tidal  displacement,  and  tube 
diameter  are  known.  This  inflection  point  will  provide  the 
information  to  deduce  the  density  and  kinematic  viscosity. 

Both  of  the  above  methods  assume  that  the  gas  under 
test  does  not  have  any  change  of  properties  along  the  tube 
length  during  the  test,  which  is  difficult  to  attain 
experimentally. 

With  the  non-dimensional  parameters  and  theoretical 
functions  derived,  the  experimental  data  can  be  evaluated. 
The  results  of  the  experiment  are  presented  in  the  next 
chapter. 


CHAPTER  5 


EXPERIMENTAL  RESULTS 
Setup  One 

The  first  experimental  setup  was  run  using  oxygen  as 
the  carrier  gas  and  nitrogen  in  air  as  the  upstream 
diffusing  gas.  These  gases  were  chosen  as  best  for  this 
experiment  because  the  densities  and  kinematic  viscosities 
did  not  vary  between  the  gases  appreciably  and  were  gases 
readily  available.  The  raw  data  will  be  presented  as  well 
as  the  computed  results.  Two  isolators  (end  chambers)  were 
available  during  these  experiments.  These  isolators 
separated  the  pump  from  the  tube,  insuring  that  no  gas 
leaked  into  or  out  of  the  pump.  The  data  was  collected 
using  the  larger  of  the  two  isolators,  due  to  the  isolator- 
tube  induced  harmonics  at  low  frequencies  for  the  smaller 
isolator. 

Achieving  steady  state  conditions  in  some  of  the  test 
runs,  sometimes  required  over  six  hours.  Because  of  this 
long  time  span,  it  was  difficult  to  determine  if  a steady 
state  condition  in  fact  was  achieved.  Some  deviations  in 
the  results  are  possible  because  of  non-steady  state 
conditions  being  accepted  as  steady  state.  At  times,  large 
concentration  surges  that  were  short  lived  were  observed 
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during  recording.  It  was  not  possible  to  determine  if  this 
was  due  to  turbulence,  a concentration  wave  passing  through 
the  tube,  poor  mixing  in  the  reservoirs,  or  any  other 
reason.  The  test  runs  that  had  exhibited  obvious  large 
concentration  surges  were  discarded. 

Repeating  these  experiments  using  other  carrier  gases 
was  not  possible  due  to  the  large  sampling  requirements  of 
the  available  oxygen  and  helium  concentration  meters.  Even 
the  small  sampling  requirements  of  the  nitrogen 
concentration  meter  affected  the  results  of  this  experiment. 
The  sampling  meter  could  not  be  connected  to  the  ports 
continuously,  as  this  affected  the  concentration  in  the 
tube.  Instead,  the  sampling  probe  had  to  be  inserted 
temporarily  to  obtain  a reading  and  removed  to  retain  the 
concentration  in  the  tube. 

Highly  repeatable  concentration  readings  were  not 
possible  for  an  unknown  reason,  therefore,  an  averaged 
reading  was  used  to  reduce  the  effects  of  the  variance  of 
the  readings.  It  was  felt  that  an  average  of  the  actual 
concentration  readings  was  superior  to  averaging  the  ratio 
of  the  individual  reading  since  the  concentration  readings 
were  not  instantaneous  readings,  but  reading  taken  a short 
time  interval  apart. 

The  pressure  transducer  was  calibrated  using  a water 
column,  obtaining  27.6  cm  of  water  for  30  divisions  on  an 
oscilloscope  on  a scale  of  0.5  volts  per  division.  This 
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results  in  a pressure  of  1.84  cm  of  water  per  volt  as  the 
calibration  of  the  pressure  transducer.  The  pressure 
measured  was  peak-to-peak  pressure,  therefore,  half  of  this 
value  was  used  as  the  amplitude  of  the  pressure  in  the 
pressure  gradient  measurement. 

The  three  tubes  used  were  each  153  cm  in  length.  The 
tubes  had  radius  measurements  of  0.587  cm,  0.777  cm  and  1.03 
cm.  The  kinematic  viscosity  used  for  calculations  was 
v = 0.16  cm  squared  per  second.  The  pressure  ports  were 
30.5  cm  apart,  61  cm  from  each  end  and  flush  to  slightly 
recessed  from  the  inner  surface  of  the  tube. 
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The  results  shown  in  Table  5.2  show  a relatively  large 
error,  low  for  smaller  Womersley  numbers  and  high  for  larger 
Womersley  numbers.  The  lower  values  are  possibly  in  error 
due  to  steady  state  conditions  not  occurring.  The  nearly 
consistent  error  for  the  higher  values  is  of  unknown  origin. 
Because  the  experimental  N values  are  above  all  possible  N 
function  values  for  these  Womersley  Numbers,  a minimum  error 
by  least  sum-of -squares  is  meaningless.  This  experimental 
test  run  sequence  was  thus  of  no  value  in  determining  an 
experimental  value  of  diffusion  coefficient.  The  test  setup 
was  unable  to  get  any  lower  than  a Womersley  Number  of  6.8, 
due  to  pressure  gradient  measurements  limitations. 

Table  5.2  Results  Setup  1 


Test 

Number 

Womersley 

Number 

NeXp 

N 

theory 

Nexp 

% Error 

Deffm p 

AXcalc 

cm 

1 

6.88 

0.00402 

0.00595 

-32.4 

1617 

50.86 

2 

8.22 

0.00463 

0.00511 

-9.4 

1508 

39.58 

3 

9.01 

0.00523 

0.00472 

10.8 

627 

22.26 

4 

17.08 

0.00277 

0.00262 

5.7 

471 

26.41 

5 

20.41 

0.00260 

0.00221 

17 . 6 

606 

26.23 

6 

11.93 

0.00415 

0.00367 

13.1 

1940 

45.64 

7 

11.93 

0.00412 

0.00367 

12 . 3 

1265 

37.00 

8 

13.77 

0.00370 

0.00321 

15.3 

187 

13.22 

9 

15.40 

0.00338 

0.00289 

17.0 

450 

19.38 

We  note,  however,  that  quite  large  values  of  Deff/DM 
where  found.  Such  large  increases  support  the  theoretical 
predicted  values. 
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Figure  5.1  Experimental  Data  Setup  1 
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Setup  Two 

In  the  second  setup,  a capillary  system  of  76  glass 
tubes  in  a close-packed  bundle  each  0.38  cm  in  diameter  were 
used.  These  tubes  emptied  into  reservoirs  on  each  end  which 
had  inlet  and  outlet  connections.  The  tidal  displacement 
was  measured  using  a spirograph  which  was  connected 
separately  from  the  actual  test  to  reduce  the 
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compressibility  effects  from  the  additional  load  on  the 
test.  The  flow  rate  of  each  connection  had  to  be  measured, 
attempting  to  minimize  the  net  flow  in  the  tubes  themselves. 
All  of  these  measurements  were  conducted  in  static 
conditions.  The  settling  time  for  this  experimental  setup 
was  much  quicker  than  the  first  setup,  being  on  the  order  of 
five  minutes.  This  was  due  to  the  smaller  diameter  tubes, 
the  larger  total  tube  area,  smaller  reservoirs,  and  larger 
feed  rates  of  the  gases.  Table  5.3  presents  the  raw  data 
that  was  collected  during  the  operation  of  the  second  setup. 
Table  5.3  Raw  Data  Setup  2 


Test 

Run 

Freq 

Hz 

ft 

ml/s 

C0 

% 

c, 

% 

AX 

cm 

10 

5 

500 

50.6 

35.3 

25.1 

11 

5 

520 

51.4 

35.3 

25.1 

12 

10 

520 

46.6 

41.4 

33 . 0 

13 

1 

520 

74 . 0 

5.0 

22.7 

14 

1 

520 

72.5 

4.9 

22.7 

15 

10 

530 

47.6 

41.8 

33 . 0 

16 

2.5 

500 

55.9 

20.0 

23 . 2 

Table  5.4  presents  the  calculated  results  with  the 
error  in  each  run.  The  values  of  Womersley  number  is  based 
on  the  kinematic  viscosity  of  the  carrier  fluid  being 
entirely  carbon  dioxide  (0.085  cm2/s)  , and  the  theoretical 
lambda  value  based  on  the  Schmidt  number  being  the  reference 
value  of  0.518,  taken  from  tables  2.1  and  2.2. 
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Table  5.4  Results  Setup  2 


Test 

Number 

Womersley 

Number 

^exp 

^theory 

^exp 

% Error 

Deffe>v 

10 

3 . 653 

0.0135 

0.0127 

6.3 

1629 

11 

3 . 653 

0.0133 

0.0127 

4.7 

1605 

12 

5.166 

0.0140 

0.0128 

9.4 

5841 

13 

1.634 

0.0027 

0.0035 

-22.9 

53 . 3 

14 

1.634 

0.0027 

0.0035 

-22.9 

53 . 3 

15 

5.166 

0.0129 

0.0128 

0.8 

5382 

16 

2 . 583 

0.0076 

0.0081 

-6.2 

392 

Figure  5.2  Results  Setup  2 
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The  large  error  in  test  runs  13  and  14  is  most  likely 
due  to  an  error  in  concentration  measurement  at  the  low 
concentration  reading.  Because  there  is  a relatively  small 
diffusional  flux  in  these  test  runs,  an  error  in  the 
concentration  reading  will  be  greatly  amplified.  Another 
likely  reason  is  because  of  the  great  concentration  change 
over  the  length  of  the  tube.  If  the  actual  diffusion 
coefficient  remains  constant,  the  Schmidt  number  will  change 
by  a factor  of  nearly  two  over  the  length  of  the  tube,  due 
to  the  viscosity  change  of  the  gas  mixture. 

An  estimate  of  the  Schmidt  number  can  be  determined  by 
a least  sum-of-squares  of  error  between  the  theoretical 
curve  and  collected  data.  The  least  error  by  this  method 
was  determined  to  be  a Schmidt  number  of  0.529.  This  value 
implies  a coefficient  of  diffusion  of  0.161  cm2/s.  This 
compares  very  favorably  to  the  tabulated  value  of  0.164 
cm2/s  at  the  same  conditions.  When  reduced  to  standard 
tabulated  conditions  of  0°C  and  one  atmosphere,  the 
determined  value  is  0.136  cm2/s  compared  to  the  tabulated 
value  of  0.139  cm2/s. 

The  error  bars  shown  in  figures  5.1  and  5.2  were 
computed  using  one-half  of  the  resolution  on  each 
measurement,  calculating  the  worst  case  condition  for  high 
and  low  respectively. 
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Pressure-Displacement 


The  pressure-displacement  equation  5.1,  namely. 


pco2Ax 


2 yfl  J1  (/-IOC) 
a JQ  i/^ia) 


(5.1) 


was  verified  using  air  in  the  tube.  The  tidal  displacement 
was  measured  using  a spirograph.  This  was  calibrated  by 
injecting  into  the  system  a measured  volume  of  gas  and 
noting  the  output  voltage  from  the  spirograph.  The  pressure 
gradient  was  determined  by  measuring  the  pressure  in  the 
reservoir  and  dividing  this  by  the  length  of  the  tube. 
Although  this  is  not  a true  measurement  of  pressure 
gradient,  the  tube  entry  losses  were  estimated  to  be 
minimal,  with  the  greatest  pressure  drop  from  the  flow  in 
the  tube.  The  pressure  transducer  was  calibrated  with  a 
water  column.  The  density  used  was  a tabulated  value  of 
1.177  Kg/m3.  The  volume  is  listed  as  divisions  noted  on  the 
oscilloscope  with  a conversion  factor  of  15.7  cm/div.  The 
pressure  gradient  is  also  listed  in  divisions  peak-to-peak 
as  noted  on  the  oscilloscope  with  a conversion  factor  of 
0.0000266  atm/(cm  div) . With  the  exception  of  the  last 
entry,  all  data  was  collected  in  the  second  setup.  The 
minimum  Womersley  number  that  was  obtainable  with  reasonable 
results  in  setup  one  was  about  6.8,  due  to  the  method  of 
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obtaining  Ax.  Ax  was  measured  by  a pressure  increase  in  a 
large  bottle. 


Table  5.5  Pressure-Displacement  Results 


Hz 

Ax 

div 

Press . 
Grad, 
div 

Womersley 

pto2Ax 

2 f 

dx 

Theory 

1 

1.45 

0.3 

1.203 

0.1308 

0.1757 

5 

1.6 

1.9 

2 . 691 

0.5698 

0.5703 

10 

2 . 1 

7.9 

3 .806 

0.7195 

0.6935 

3 . 5 

* 

36 

6.793 

0.8050 

0.8126 

Figure  5.3  Experimental  Data  Pressure-Displacement 
Relationship 
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CHAPTER  6 


CONCLUSIONS 
Concluding  Remarks 

This  experiment  proved  that  the  Molecular  Diffusion 
Coefficient  in  gases  can  be  determined  by  means  of  the 
enhanced  dispersion  in  oscillatory  flow.  Ideally,  data  for 
a large  range  of  Womersley  numbers  can  be  collected  and  this 
will  provide  a very  precise  determination  of  the  Schmidt 
number  from  which  the  molecular  diffusion  coefficient  then 
is  found.  This  ensures  that  the  one  correct  lambda,  or  N 
curve  as  a function  of  Schmidt  number  is  selected  and 
reduces  the  statistical  standard  deviation  from  a large 
sampling.  Once  this  Schmidt  number  is  found,  the  molecular 
diffusion  coefficient  can  be  determined  by  dividing  the 
carrier  fluid's  kinematic  viscosity  by  the  Schmidt  number. 

In  the  discussion  here,  lambda  will  be  used,  but  the  same 
conclusions  also  apply  to  the  N function. 

If  only  a small  range,  or  small  sampling  of  Womersley 
number  can  be  obtained,  from  either  eguipment  limitations  or 
any  other  causes,  the  range  of  collected  data  will  be  most 
sensitive  for  Womersley  number  near  the  tuning  point.  That 

d\ 

is,  the  most  sensitive  data  occurs  near  = 0,  which 
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corresponds  to  the  peak  of  the  lambda  function.  This  is  so 
critical  because  it  will  allow  a check  of  the  measured  value 
of  lambda  itself  which  is  not  possible  on  other  limited 
ranges  of  the  lambda  curve.  If  sufficient  data  is  collected 
to  specify  the  actual  peak  of  this  curve,  the  appropriate 
Schmidt  number  can  be  determined  directly. 

The  second  best  limited  range  for  gases  is  in  the  range 

dX 

of  Womersley  numbers  below  the  peak,  where  -^>0.  This 

selection  is  desirable  because,  as  a function  of  Schmidt 
number,  the  lambda  function  near  the  desired  Schmidt  number 
is  always  monotomically  increasing.  This  allows  for  far 
better  discrimination  of  the  Schmidt  number  using  the  least 
sum-of-squares  of  error  procedure. 

The  least  desirable  location  for  data  points  is  above 

the  peak  of  the  lambda  curve,  or  <0.  In  this  location, 

the  function  is  not  necessarily  monotomic  as  a function  of 
Schmidt  Number. 

The  greatest  problem  in  determining  the  diffusion 
coefficient  with  this  method  is  in  determining  the  tidal 
displacement.  Because  this  term  is  squared  in  the 
experimental  determination  of  lambda,  any  error  is  greatly 
amplified.  In  practice,  the  determination  of  N should  be 
more  precise  than  lambda,  due  to  the  generally  more  precise 
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measurement  of  the  magnitude  of  pressure  gradient  instead  of 
tidal  displacement. 

Although  this  procedure  appears  to  work  well  with 
gases,  in  liquids  it  may  be  far  less  desirable.  Liquids 
generally  have  Schmidt  numbers  on  the  order  of  1000,  whereas 
in  most  gases  it  is  near  0.7.  With  the  Schmidt  number  this 
large,  the  lambda  curve  peaks  well  below  a Womersley  number 
of  one.  In  this  range,  there  is  very  little  distinction 
between  different  lambda  curves  as  a function  of  Schmidt 
number,  therefore,  a precise  determination  of  Schmidt  number 
becomes  much  more  difficult. 

It  is  interesting  that  the  presented  method  for 
determining  DM  is  found  to  be  quite  accurate  for  gases  and 
thus  offers  an  alternative  to  the  standard  Loschmidt 
approach  for  determining  this  important  constant 
experimentally . 


Future  Work 


Future  work  should  involve  a study  of  the  effect  of 
harmonics  on  the  observed  enhanced  dispersion  considered 
here.  The  collected  data  with  setup  one,  suspects  that 
harmonics  may  increase  the  effective  diffusion  in  the  range 
of  Womersley  Number  greater  than  the  tuning  point  of  the  X 


curve . 
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Some  consideration  should  also  be  given  to  the  use  of 
transient  changes  in  concentration  at  a specific  location 
within  the  confining  tubes.  Such  an  approach  would  allow 
not  only  a more  rapid  determination  of  DM  but  also  would 
allow  extension  of  the  present  method  to  high  Schmidt  number 
fluids  such  as  liquids. 


APPENDIX 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c C 

C CALCULATES  LAMBDA  GIVEN  THE  WOMERSLEY  AND  SCHMIDT  NUMBERSC 

C AND  C 

C CALCULATES  N GIVEN  THE  WOMERSLEY  AND  SCHMIDT  NUMBER  C 

C c 

C WOMERSLEY  NUMBER  = RADIUS  OF  TUBE  TIMES  SQUARE  ROOT  OF  C 
C QUANTITY  ANGULAR  FREQUENCY  OVER  C 

C DYNAMIC  VISCOSITY  C 

C C 

C SCHMIDT  NUMBER  = DYNAMIC  VISCOSITY  OF  CARRIER  GAS  C 

C OVER  MOLECULAR  DIFFUSION  COEFFICIENT  C 

C OF  TRACE  GAS  INTO  CARRIER  GAS  C 

C C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
DOUBLE  PRECISION  FAC(0:80),  A,  S,  LAM,  N 
INTEGER  I 
COMMON  FAC 

WRITE ( * , * ) ' Enter  a value  for  Womersley  Number  (A) ' 

READ ( * , * ) A 
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80 


WRITE (*,*) 'Enter  a value  for  Schmidt  Number  (S) ' 

READ ( * , * ) S 
FAC ( 0 ) = 1D0 
DO  10  I = 1,  80 
FAC ( I ) = FAC(I-l)  * I 
10  CONTINUE 

IF  (SQRT(S)  * A .GT.  85D0  .OR.  A .GT.  85D0)  THEN 
WRITE (*, *) 

WRITE (*,*) 'POSSIBLE  ERROR  DUE  TO  EXCEEDING  THE  MAXIMUM 
#ALLOWABLE  FOR  THIS  APPROXIMATION' 

WRITE(*,*) 'THE  MAXIMUM  FOR  A,  OR  SQRT(S)  * A IS  85' 

ENDIF 

WRITE (*, *) 

WRITE ( * , * ) ' (Womersley)  (Schmidt) ' 

WRITE ( * , * ) ' A S LAMBDA ( A , S ) 

# N (A, S) ' 

WRITE(*, 101) A,  S,  LAM (A,  S) , N (A,  S) 

101  FORMAT ( 3X , F6.3,  5X,  F6 . 3 , 6X,  F18.16,  2X,  F18.16) 


END 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c c 

C SUBROUTINE  TO  CALCULATE  LAMBDA  AS  A FUNCTION  OF  BESSEL  C 
CFUNCTIONS  WITH  PARAMETERS  OF  WOMERSLEY  AND  SCHMIDT  NUMBERSC 
C C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

DOUBLE  PRECISION  FUNCTION  LAM (A,  S) 

DOUBLE  PRECISION  A,  S 

DOUBLE  COMPLEX  SRNI,  J,  JO,  Jl,  JOS,  J1S 
SRNI  = CDSQRT ( (0D0,-1D0)  ) 

JO  = J(0,  SRNI  * A) 

Jl  = J(l,  SRNI  * A) 

JOS  = J(0,  SRNI  * DSQRT(S)  * A) 

J1S  = J(l,  SRNI  * DSQRT(S)  * A) 

LAM  = DBLE ( ( -DSQRT(S)  * Jl  * JOS  / (JO  * J1S)  + S)  * 
#CDSQRT(  ( 0D0 , IDO ) ) * DCONJG(Jl)  / DCONJG(JO))  / 

#(4D0  * A * (IDO  - S **  2)  * ( CDABS ( IDO  - 2DO  * Jl  / 

# ( SRNI  * A * JO) ) **  2) ) 

RETURN 


END 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c C 

C SUBROUTINE  TO  CALCULATE  N AS  A FUNCTION  OF  BESSEL  C 

CFUNCTIONS  WITH  PARAMETERS  OF  WOMERSLEY  AND  SCHMIDT  NUMBERSC 
c C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

DOUBLE  PRECISION  FUNCTION  N (A,  S) 

DOUBLE  PRECISION  A,  S 

DOUBLE  COMPLEX  SRNI,  J,  JO,  Jl,  JOS,  J1S 
SRNI  = CDSQRT ( (0D0,  -IDO)  ) 

JO  = J(0,  SRNI  * A) 

Jl  = J(l,  SRNI  * A) 

JOS  = J(0,  SRNI  * DSQRT(S)  * A) 

J1S  = J(l,  SRNI  * DSQRT(S)  * A) 

N = DBLE ( ( -DSQRT(S)  * Jl  * JOS  / (JO  * J1S)  + S)  * 

#CDSQRT ( ( 0D0 , IDO ) ) * DCONJG(Jl)  / DCONJG(JO))  / 

# ( 4D0  * A * (IDO  - S **  2)  ) 

RETURN 


END 


83 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c C 
C SERIES  REPRESENTATION  OF  THE  BESSEL  FUNCTION  C 
C C 
C V = ORDER  OF  FUNCTION  C 
C X = POINT  TO  BE  EVALUATED  C 
C C 
C VALID  FOR  ABS(Y)  < 85  WHEN  X = Y * SQRT(-I)  C 
C C 
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


DOUBLE  COMPLEX  FUNCTION  (V,  X) 

DOUBLE  COMPLEX  X,  J1 
DOUBLE  PRECISION  FAC (0:80) 

INTEGER  V,  I 
COMMON  FAC 
J1  = ( 0D0 , 0D0 ) 

DO  10  I = 0,79 

J1  = J1  + ( -IDO)  **  I * (X  / 2D0 ) **  (2  * I + V)  / 
# ( FAC ( I ) * FAC (I  + V) ) 

10  CONTINUE 
J = J1 
RETURN 


END 
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